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Abstract 

We calculate the initial non-equilibrium conditions from perturbative QCD 
(pQCD) within Glauber multiple scattering theory for y/s = 200 AGeV and 
\fs = 5.5 ATeV. At the soon available collider energies one will particularly test the 
small x region of the parton distributions entering the cross sections. Therefore 
shadowing effects, previously more or less unimportant, will lead to new effects 
on variables such as particle multiplicities dN/dy, transverse energy production 
dEx/dy, and the initial temperature Tj, In this paper we will have a closer look 
on the effects of shadowing by employing different parametrizations for the shad- 
owing effect for valence quarks, sea quarks and gluons. Since the cross sections 
at midrapidity are dominated by processes involving gluons the amount of their 
depletion is particularly important. We will therefore have a closer look on the 
results for dN/dy, dEx/dy, and Tj by using two different gluon shadowing ra- 
tios, differing strongly in size. As a matter of fact, the calculated quantities differ 
significantly. 
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1. Introduction 

One of the challenging goals of heavy ion physics is the detection of the 
quark-gluon plasma, a state in which the partons are able to move freely 
within a distance larger than the typical confinement scale r co „/. ~ 1/Aqcd ~ 
1/0.2 GeV ~ 1 fm. The build-up of this state should happen early in a 
heavy ion reaction when the two streams of initially cold nuclear matter pass 
through each other. Thereby first virtual partons are transformed to real ones 
and later on in the expansion phase the fragmentation of the partons into col- 
orless hadrons takes place. When separating pQCD from non-perturbative 
effects at some semi-hard scale Pq = 2 GeV the respective time scale of 
perturbative processes is thus of the order r ~ 1/po ~ 0.1 fm/c which ap- 
proximately coincides with the lower bound of the initial formation time of 
the plasma in a local cell [|IJ . Therefore all further evolution of the system is 
significantly influenced by the initial conditions of pQCD since macroscopic 
parameters, as e.g. the initial temperature Tj, directly enter into hydrody- 
namical calculations. 

We here will focus on the very early phase of an ultrarelativistic heavy ion 
collision and use pQCD above the semi-hard scale p s h. = Po = 2 GeV. 

In a typical high energy pp or pp event one measures distinct hadronic 
jets with a transverse momenta of several GeV (pt > 5 GeV) 0. In con- 
trast to the experimental very clean situation of hadronis jets at large p? one 
encounters the problem of detectability of low transverse momentum jets in 
heavy ion collisions. These so-called minijets contribute significantly to the 
transverse energy produced in AB collisions due to their large multiplicity 
[[J. The major part of these set-free partons are gluons that strongly domi- 
nate the processes as their number is much larger for the relevant momentum 
fractions. In turn the shadowing effects are expected to be much larger for 
gluons than for the quark sea 0]. Therefore the relative contribution of the 
gluons should decrease but still dominate the cross sections. The shadowing 
of the gluons has the peculiarity of not being known exactly due to the neu- 
trality of the mediators of the strong interaction which makes it impossible 
to access Rg(x, Q 2 ) directly in a deep inelastic e + A event. Therefore we 



2 



will here investigate two possible parametrizations of the shadowing ratio 
Rg = xG A I A ■ xG N for gluons as will be described below in detail. 

2. Minijets 

As outlined above, we will here investigate the effects of shadowing on the 
minijet production cross sections. The production of a parton / = g,q,q can 
in leading order be described as |§ 



dy 



J dp 2 T dy 2 J2 x ifi( x i>Q 2 ) x 2 fj(x 2 ,Q 2 

ij,kl 

' da l ^ kl r A . p da^ kl _ . 
df k — -j — (t,u) + bfi — Tl — (u,t) 



dt 



dt 



1 + 5, 



ki 



(1) 



The factor 1/(1 + Ski) enters due to the symmetry of processes with two 
identical partons in the final state. The exchange term da(t, u) <-> da(u, t) 
accounts for the possible symmetries of e.g. having a quark from nucleon i 
and a gluon from nucleon j and vice versa, i.e. it handles the interchange of 
two of the propagators in the scattering process. The possible combinations 
of initial states are 



*J = 99, 91, Q9, 9Q, 19, QQ, QQ, QQ, QQ 



(2) 



The momentum fractions of the partons in the initial state are 
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The integration regions are 
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M<ln (f + 
V 2 Po 



Ap 2 



The mandelstam variables are defined as 



(5) 



X\ ■ X2 ■ S, t 



-Pt 



1 + e 



(2/2 -y) 



u = —p T 
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For the parton distributions entering the handbag graph we choose the GRV 
LO set H for RHIC. Since at LHC one probes smaller momentum fractions 
we there use the newer CTEQ4L parametrization |J with Nf = 4 and Q = 
Pt- The normalization is done so that one has two outgoing partons in one 
collision, i.e. 

dy-J- = 2 °Lrd (7) 



dy 

In the calculations the boundaries for the calculations are either over the 

whole rapidity range or \y\ < 0.5 for the central rapidity region. 

To account for the higher order contributions at RHIC we choose a fixed K 

factor of K=2.5 from comparison with experiment as discussed in 0, [|. In 

the range 5.5 GeV < px < 25 GeV a factor K=2.5 is needed to describe 

the UAl data, and in the range 30 GeV < px < 50 GeV a factor of K=1.6 

is needed. However the cross section has dropped so much at these large 

transverse momenta that we keep K=2.5 fixed for all pt- For LHC energies 

the mean px tends to be larger; so we choose K=1.5 for this case. 

By applying Glauber theory we calculate the mean number of events per unit 

of rapidity: 

™L = TAA w*$=i ( 8) 

ay ay 

where the nuclear overlap function 7^(6) for central events is given by 
Taa{0) ~ A 2 /ttR a . For the nuclei in our calculation this gives Ta u a u {Q) = 
29 jmb and TpbPb(0) = 32 /mb. Again it should be emphasized that dN^ /dy 
gives the number of collisions and that the total number of partons is as 
twice as high in a 2 — ► 2 process. The necessary volume, needed to derive 
the densities from the absolute numbers, is calculated as 

V i = 7iR 2 A Ay/p , Ra = A 1 ^ x 1.1 fm (9) 

Therefore we get V { (Au + Au) = 12.9 fm 3 and Vi(Pb + Pb) = 13.4 fm 3 . 
For the energy density at midrapidity we need the first Et moment: 

da f 
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Figure 1: Rp 2 vs. R G at Q 2 = 4 GeV 2 for m Pb 



J dp 2 T dy dy 2 J2 x ifi( x iiQ 2 ) x 2fj(x 2 ,Q 2 



ij,kl 



X 



dt 



da ij ^ kl 

-(t,u) + 6 fl -j — (u,t) 

dt 



1 + 6, 



p T e(y) (10) 



Here the acceptance function e(y) is e(y) = 1 for \y\ < 0.5 and e(y) = 
otherwise. 

3. Nuclear Shadowing 

In heavy ion collisions one has to account for an effect that does not appear 
for processes involving two nucleons only: nuclear shadowing. In the lab 
frame the deep inelastic scattering at small Bjorken x (x <ti 0.1) proceeds 
via the vector mesons as described in the vector meson dominance model 
(VMD) where the handbag graph contribution becomes small. In VMD 
the interaction of the virtual photon with a nucleon or nucleus is described 
as a two step process: the photon fluctuation into a qq pair (the p,uj,(p 
mesons at small Q 2 ) within the coherence time l c and a subsequent strong 
interaction with the target ||. The coherence time arises in this picture from 
the longitudinal momentum shift between the photon and the fluctuation: 
l c ~ 1/AA; 2 where /S.k z = k] — k h z . The cross section is: 



(x(7*A0 



dz Id r 



iz,r) 



qqN 



(11) 
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Figure 2: R F , 2 vs. R G at Q 2 = A GeV 2 for 40 C a. 



where the Sudakov variable z gives the momentum fraction carried by the 
quark or the antiquark. The interaction of the fluctuation with the nucleon 
can be described in the color transparency model as M 



7T 



e qm ='-r 2 a s {Q> 2 )x>g{x\Q' 2 ) 



(12) 



where x' = M 2 g/(2mi/), r is the transverse separation of the pair and Q' 2 = 
4/r 2 . For the interaction of the fluctuation with a nucleus one makes use of 



Glauber-Gribov multiple scattering theory [10] where the fluctuation inter- 
acts coherently with more than one nucleon in the nucleus when the coherence 
length exceeds the mean separation between two nucleons: 



<r qqA = J d 2 b (l - 



-<y q q N T A {b)/2 



(13) 



When expanding for large nuclei and taking the dominating double scatter- 
ing term only one finds 



a hA = Aa 



l,N 



1 - A 1 '^ + 



(14) 



with a = 1.1 fm. 

Figures [I] and show the results for 207 Pb and 40 C a (for further details see 
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A very different scenario is employed in parton fusion models. Here the 
process of parton parton fusion in nuclei can be understood as an overlapping 
of quarks and gluons that yields a reduction of number densities at small 
x and a creation of antishadowing for momentum conservation at larger x 
The onset of this fusion process can be estimated to start at values 
of the momentum fraction where the longitudinal wavelength (1/xP) of a 
parton exceeds the size of a nucleon (or the inter-nucleon distance) inside 
the Lorentz contracted nucleus: 1/xP ~ 2R n M n / P, corresponding to a value 
x ~ 0.1. Originally the idea of parton fusion was proposed in [0 and later 



proven in to appear when the total transverse size 1/Q of the partons 
in a nucleon becomes larger than the proton radius to yield a transverse 
overlapping within a unit of rapidity, xG(x) > Q 2 R 2 . The usual gluon 
distribution in the nucleon on the light cone in light-cone gauge (n ■ A = 
A + = 0) is given by 

x G{x) = -{n-f J^-(P \F + »(0)F + ^(Xn) \ P) (15) 

The recombination is then described as the fusion of two gluon ladders into a 
single vertex. One finally arives at a modified Altarelli-Parisi equation where 
the fusion correction enters as a twist four light cone correlator. Typically 
the fusion correction in the free nucleon turns out to be significant only for 
unusually small values of x or Q 2 . As shown in [14| the situation changes 
dramatically in heavy nuclei. Here the strength of the fusion for ladders com- 
ing from independent constituents increases and is of the same order as the 
fusion from non-independent constituents. Therefore, parton recombination 
is strongly increased in heavy nuclei of A ~ 200. 

Unfortunately the different models do not give the same results for the 
ratio Rg{x,Q 2 ). We will therefore use two versions of parametrizations to 
investigate the effects of shadowing on the relevant variables. On the one 
hand we use a Q 2 dependent version of Eskola, Kolhinen, Salgado, and Ru- 
uskanen, often referred to as "'98 shadowing" (see figure |3|), that tries to 
avoid any model dependence by using sum rules for baryon number and 



momentum [15| and on the other hand we use a modified version of a Q in- 
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Figure 3: Shadowed parton distributions as parametrized by Eskola et al 
T5[ . Note that the gluon shadowing appears to be weaker than quark 



in 



shadowing and that the onset happens for smaller x. 



dependent parametrization (see figure |j) given in 1 16] which employs a much 
stronger gluon shadowing in accordance with the results of 0]. Especially 
for RHIC, where the lower bound for the momentum fraction at midrapidity 
for p T = p = 2 GeV is given by x = 2p T / ^fs = 0.02, the onset of the gluon 
shadowing, i.e. the transition region between shadowing and antishadowing, 
is of great importance. 

In |]T3| the onset of gluon shadowing (Rq = 1) is chosen at x ~ 0.029 for 
Q = 2 GeV motivated by the results found in [T7J where the connection be- 
tween the gluon distribution and the Q 2 dependence of F 2 via the DGLAP 
equations was employed: 



<91n Q 2 



J2e 2 ixG(2x,Q 2 ) 



(16) 



By using the NMC data |18| on deep inelastic scattering on a combina- 
tion of Sn and C targets the ratio G Sn (x)/G c (x) was derived in the range 
0.011 < x < 0.18. The cross over point can, despite the large errorbars, be 
guessed to be x ~ 0.03. However one should add here that the situation 
for Rq 3 = xG Pb (x) I xG N (x) can look rather different. Since this question of 
the onset of gluon shadowing is not yet settled we chose the same onset for 
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1.2 




Figure 4: Our variation of the parametrization given in [15 with much 
stronger gluon shadowing as found in the VMD calculation in ||. The 
stronger shadowing is also motivated by the fact that we calculate central 
collisions where the shadowing effect is stronger than for 6-averaged collisions. 



quark and gluon shadowing in our modified parametrization to investigate 
the relevance of this point. We fixed Rp 2 = Rg = 1 at x ~ 0.07. In VMD as 
well as in parton fusion models the onset is treated on an equal footing: for 
the coherent scattering processes in VMD it should make no difference (at 
least for the onset) whether a qq or a gg pair scatters from more than one 
nucleon at l c > r^^. In the parton fusion model one treats the leaking out 
of the partons equally for sea quarks and for gluons since for both sea quarks 
and for gluons one has a spatial extent of 1/ xP in the longitudinal direction 
and therefore the onset for Rg and Rf 2 is essentially the same in this model. 

4. Results 

In the following we will give the results for the different parton species 
/ = g, q, q at RHIC and LHC including the different shadowing parametriza- 
tions or none shadowing, respectively. The results for the number of partons 
J dN-f jdy can easily be derived from / dy da? jdy by the relation dN? /dy = 
TAA{Q)d<jf jdy. All results include a K-factor of K=2.5 for RHIC and K=1.5 
for LHC. On the one hand we give the results for the whole y-range and 
on the other hand we give the result for the central rapidity region which is 
of special interest, not only from the experimental setup point of view but 
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also since it is the region where highest parton densities and the strongest 
shadowing effects are expected. 

Let us start by giving the results without shadowing corrections for RHIC. 
The first three tables give the unshadowed multiplicities integrated over the 
whole rapidity range and over the central region, respectively. Tables 4 
through 6 give the first E? moments for the respective parton species. The 
rapidity distributions for the cross sections are depicted in figure ||. 



Table 1: / dy dW/dy for yfs = 200 AGeV 



range of y 


99 -> 99 


gq-*gq + gq-+ gq 


TOTAL 


all y 


920.8 


384.3 


1305.1 


\y\<0.5 


192.9 


90.7 


283.6 



Table 2: J dy dN q /dy for y/j = 200 AGeV 



range of y 


gq -> gq 


qq -> qq 


gg -> qq 


qq^qq 


TOTAL 


all y 


310.3 


57.3 


6.5 


22.8 


396.9 


\y\ < 0.5 


21.0 


7.4 


1.5 


2.3 


32.2 


Table 3: / dy dW/dy for = 200 AGeV 


range of y 


9<i-> gq 


qq^qq 


gg -» qq 


qq^qq 


TOTAL 


all y 


74.2 


22.8 


6.5 


2.2 


105.7 


\y\ < 0.5 


12.5 


5.3 


1.5 


0.4 


19.7 
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Figure 5: Unshadowed rapidity distributions of gluons, quarks, and anti- 
quarks. 



11 



Table k.cjv (E T ) [mb GeV] 



range of y 


55 -»• 55 


gq ^ gq + gq ^ gq 


TOTAL 




y 


< 0.5 


18.02 


8.72 


26.74 



Table h:a q (E T ) [mb GeV] 



range of y 


gq -> #g 


gg -> gg 


55 -»• 99 


gg -> gg 


TOTAL 




y 


< 0.5 


2.065 


0.786 


0.1398 


0.22 


3.2 



Table 6:<r g ~ (E T ) [mb GeV] 



range of y 


gq -> #g 


gg -> gg 


55 -»■ 99 


gg -> gg 


TOTAL 




y 


< 0.5 


1.206 


0.51 


0.139 


0.004 


1.896 
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Figure 6: Rapidity distributions of gluons, quarks, and antiquarks with our 
modified shadowing parametrization shown in figure 
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For the strong gluon shadowing shown in figure |] one finds the following 
multiplicities for the different parton species (the rapidity distributions are 
shown in figure |5|) : 



Table 7: J dy dN^/dy for y/s = 200 AGeV 



range of y 


99 -> 99 


gq -> gq + gq -> #9 


TOTAL 


all y 


581.5 


249.4 


830.9 


y\ < 0.5 


122.6 


60.2 


182.8 



Table 8: J rfy oWyrfy /or = 200 AgeF 



range of y 


#9 -> £9 


gg -> gg 


gg -> 


gg -> gg 


TOTAL 


all y 


196.5 


48.6 


3.6 


18.1 


266.8 


< 0.5 


15.9 


5.9 


0.9 


1.7 


24.4 


Table 9: J rfy diV«/dj/ /or ^ = 200 AGeV 


range of y 


gq^ gq 


qq-^qq 


-> 99 


qq-^qq 


TOTAL 


all y 


52.9 


18.1 


3.6 


1.8 


76.4 


|S/| < 0.5 


9.4 


4.1 


0.9 


0.3 


14.7 



The first moments for the reactions including our modified strong gluon 
shadowing are given by: 

Table 10: (E T ) [mb GeV] 



range of y 


gg -> 


gq ^ gq + gq ^ gq 


TOTAL 


\y\<0.5 


11.87 


5.93 


17.8 
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Table 11: a q (E T ) 



mb GeV] 



range of y 


gq -> gq 


qq -> gg 


99 


gg -> gg 


TOTAL 


\y\ < 0.5 


0.64 


0.25 


0.04 


0.072 


1.002 



Table 12: of (E T ) [mb GeV] 



range of t/ 


gq -> #9 


gg -> gg 


-> 99 


gg -> gg 


TOTAL 




2/ 


< 0.5 


0.37 


0.013 


0.037 


0.165 


0.585 



We also calculated the multiplicities and first £r moments by employing the 
newest available shadowing parametrization of Eskola et al of ref. JT5] shown 
in figure [3|. As emphasized above one should note that the shadowing of 
gluons in this parametrization is smaller than the quark shadowing since it 
was tried to stay away from any model dependence and just stick to sum 
rules expressing the momentum and baryon number conservation but still 
assuming that at small x (x ~ 10~ 4 ) the gluon ratio should coincide with the 
sea quark ratio. By employing this version we find the results listet in the 
following tables and shown in figure |7| : 



Table 13: J dy dN^/dy for y/s = 200 AGeV 



range of y 


gg -> gg 


gq -> gq + gq -> gq 


TOTAL 


all y 


969.5 


350.2 


1319.7 


y < 0.5 


201.8 


81.9 


283.7 



Table 14: 



dy dN q /dy for y/s = 200 AGeV 



range of y 


gq -> gq 


gg -> gg 


gg -> 99 


gg -> gg 


TOTAL 


all y 


282.8 


50.8 


65.3 


18.9 


417.8 


y\ < 0.5 


19.6 


6.3 


1.5 


1.7 


29.14 
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Figure 7: Rapidity distributions of gluons, quarks, and antiquarks with the 
'98 version of the shadowing parametrization shown in figure |3|. 



Table 15: J dy dN^/dy for y/s = 200 AGeV 



range of y 


gq -* gq 


qq^qq 


99 -> (R 


qq-^qq 


TOTAL 


all y 


66.8 


18.7 


65.3 


16.7 


167.5 


\y\<0.5 


10.9 


4.1 


1.5 


0.3 


16.9 
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Figure 8: Comparison of the rapidity distributions with strong gluon shadow- 
ing (left figure) and weak shadowing (right figure) to unshadowed distribution 
for RHIC (see text). 

In figure || we directly compared the strong gluon shadowed distributions 
(left figure) with the unshadowed one. The same was done for the compar- 
ison of the Q 2 dependent '98 shadowing version [[OJ with the unshadowed 
one (right figure). The solid lines give the total contribution, the dotted 
ones the contribution from the gg subprocess and the dashed lines give the 
gq + gq contribution. The thick lines denote the unshadowed distributions 
and the thin ones the two shadowed ones. Note that due to the onset of 
gluon shadowing in the '98 version at such small values of x one even gets an 
enhancement for the gg — > gg subprocess at RHIC. We also calculated the 
Pt distribution without and with the two shadowing versions at midrapid- 
ity (figure |9|). Unlike the strong shadowing case the cross over point of the 
curves already happens at pr ~ 2.5 GeV for the '98 gluon shadowing version 
which immediately explains the enhancement in the rapidity distribution. 
For the first Ex moment of the transverse energy we find with the shadowing 
parametrization of Eskola et al 



Table 16: a 9 (E T ) [mb GeV] 



range of y 


99 -> 99 


gq -> gq + gq -> gq 


TOTAL 


\y\<0.5 


19.2 


8.05 


27.25 
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Figure 9: px distributions for the two shadowing parametrizations at midra- 
pidity. 



Table 17: a q (E T ) 



mb GeV] 



range of y 


gq -> gq 


qq -> gg 


-> <?<? 


qq -> gg 


TOTAL 




y 


< 0.5 


2.03 


0.67 


0.15 


0.018 


2.87 



Table 18: (E T ) [mb GeV] 



range of y 


#g -> #g 


gg -> gg 


-> 


gg -> gg 


TOTAL 




y 


< 0.5 


1.123 


0.032 


0.148 


0.423 


1.726 



From the results above we can calculate the total transverse energy Ej> — 
a < Et > Taa(0) carried by the partons, the number and energy densi- 
ties rif and £/, and also derive the initial temperature Tj if we assume the 
behavior of an ideal gas of partons. To do so we need the initial volume. 
With R A = A l/Z ■ 1.1 fm, T AuAu {0) = 29/mb, and R Au = 6.4 fm we find 
Vi = 7iR\Ayr = 12.9 fm 3 . 

Therefore at RHIC without any shadowing and with K=2.5 we have at midra- 
pidity a total number of 284 gluons, 32 quarks, and 20 antiquarks. These 
carry a transverse energy of 774 GeV (gluons), 93 GeV (quarks), and 55 
GeV (antiquarks). 

It is then straight forward to derive the number densities by dividing by the 
initial volume to yield: n g = 22 fm~ 3 , n q = 2.5 fm -3 , riq = 1.5 fm~ 3 . The 
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energy densities can be derived in an analogous way to give: 
e s = 60 GeV/fm 3 , £ q = 7.2 GeV/fm 3 , and = 4.3 GeV/fm 3 . If we 
assume total equilibrium we can derive the initial temperature from these 
numbers as 

£tdeal = 167r 2 A T 4 (17) 

At this point some comments are appropriate: one could wonder whether 
the system can be in equilibrium since one has only hard 2^2 parton scat- 
terings in this Glauber approach. Also one often assumes global equilibrium 
to be established after, say lfm/c. Now here we are mainly interested in 
local equilibrium as it is required for example for hydrodynamical calulca- 
tions. The equilibration of partons in a local cell happens to be much faster 
for the following reasons. The high Q 2 hard scatterings among the partons 
are absolutely unimportant for the equipartition of longitudinal and trans- 
verse degrees of freedom. It are the soft interactions that are responsible 
for this feature and there is a huge resource of soft partons available in the 
nucleons, even when assuming the parton distributions to be shadowed in 
heavy nuclei. The link to the short equilibration time is the fact that even 
though the nucleus is Lorentz contracted to L/cosh y, the partons obey the 
uncertainty principle and are therefore smeared out to distances 1/xP in the 
infinite momentum frame and so the major part of the partons is outside 
the Lorentz contracted disk. Based on some basic priciples and by using the 
Fokker-Planck equation [|IJ |19| the time it takes to establish local equilibrium 
in a cell was estimated to have a lower bound of r ~ 0.15 fm/c. As noted 
above we introduced a lower momentum cut-off p = 2 GeV corresponding 
to a proper time of about 0.1 fm/c. So therefore we may not be far from 
local equilibration and the calculation on the initial temperature from the 
initial energy densitiy could be rather justified. 

For the temperature we take into account only the gluons due to their 
large multiplicity and energy density that dominates the respective values 
for the quarks. We then find T ; = 549.52 MeV for RHIC. If we neglect all 
higher orders, i.e. take a K-factor of K=l (which of course is wrong, but it 
is instructive to see the impact on Tj), we get T Z K=1 = 437 MeV. 
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The same quantities were then calculated for the two different shadowing 
scenarios. For the calculations employing the strong gluon shadowing we 
found that there are 183 gluons, 25 quarks, and 15 antiquarks carry- 
ing transverse energies of 516 GeV (gluons), 29 GeV (quarks), and 17 
GeV (antiquarks). The resulting number and energy densities are found 
to be n g = 14 fm 3 , n q = 1.9 fm 3 , n q = 1.2 fin -3 , e s = 40 GeV/fm 3 , 
e q = 2.3 GeV/fm 3 , and e q = 1.3 GeV/fm 3 . When we calculate the ini- 
tial temperature for an ideal parton gas from these numbers we find that the 
initial temperature decreases due to the reduced number and energy densi- 
ties having their origin in the shadowing of the parton distributions. We find 
Ti )S had = 496.5 MeV for a K factor of 2.5 and when neglecting all higher 
order contributions we derive T^ s ^ d = 394.9 MeV. So what we can learn 
here ist the following: due to the reduced number of partons involved in 
the hard processes a reduction in the number densities and therefore in the 
energy densities entering the formula for the temperature of a thermalized 
parton gas results. One should note that the onset of shadowing in our 
modified shadowing parametrization was chosen same for quarks and gluons 
in accordance with the onset of coherent scattering of a quark antiquark or 
gluon gluon pair, respectively off a nucleus. Now in the second shadowing 
parametrization we employed \TE\ one finds that the onset of shadowing for 
gluons starts at smaller momentum fractions from xG Sn (x)/xG (x) data. 
With a momentum cut-off po = 2 GeV the momentum fractions involved 
in processes at midrapidity are bound from below at x — 0.02. Therefore 
one is right on the edge of the onset of shadowing of the parametrizations 
and one should expect the very interesting case that one is on the edge 
to the antishadowing region for gluons in the '98 parametrization of Es- 
kola et al but not so for the parametrization employing the strong gluon 
shadowing. This behavior is immediately reflected in the number and en- 
ergy densities. We found that for this specific shadowing parametrization 
one has 284 gluons, 29 quarks, and 17 antiquarks carrying transverse 
energies of 790 GeV (gluons), 83 GeV (quarks), and 50 GeV (anti- 
quarks). We found the following densities: n g = 22 fm 3 , n q = 2.2 fm -3 , 
n q = 1.3 fm 3 , e g = 61.2 GeV/fm 3 , e q = 6.43 GeV/fm 3 , and 
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£q = 3.88 GeV/fm 
T 

i,shad 



3 These numbers result in an initial temperature of 
439.2 MeV } respectively. 



552.3 MeV and T^ d 



We also went through the same program to investigate the impact of the 
different shadowing parametrizations at the higher LHC energy of \/s = 5.5 
TeV. We here used the newer parton distributions of CTEQ4L since the in- 
volved momentum fractions are so small that any new information at small 
x are valuable. When comparing GRV '94 and CTEQ4L one finds a dif- 
ference of about a factor of two at x ~ 10~ 5 . At LHC energies the effect 
of shadowing should be much more relevant than at RHIC due to the re- 
gion of smaller x that gets probed. Because of the strong dominance of the 
gluon component in the nucleon we restricted ourself to the calculation of 
a 9 , N 9 , and therefore on the transverse energy and temperature produced 
by the final state gluons only. Let us first begin with the unshadowed results. 



Table 19: J dy dN 9 /dy for s/s = 5.5 ATeV 



range of y 


99 -> 99 


gq->gq + gq-+ gq 


TOTAL 


all y 


36822.7 


6006.1 


42828.8 


\y\<0.5 


4137.6 


707.2 


4844.8 


Table 20: a 9 (E T ) [mb GeV] 


range of y 


99 -> 99 


gq -> gq + gq -> gq 


TOTAL 


y\ < 0.5 


438.09 


74.92 


513.01 



The rapidity distributions for unshadowed and shadowed gluons at LHC is 



depicted in figure 10. 
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Figure 10: Rapidity distributions of unshadowed (upper figure) and shadowed 
gluons (lower two figures) at LHC. The figure in the middle was derived by 
employing the strong gluon shadowing, whereas the bottom figure employed 
the Q 2 dependent '98 version. Note the change in shape when the strong 
gluon shadowing is employed. 
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For the strong gluon shadowing we find the following results 



Table 21: [dy dW/dy for yg = 5.5 ATeV 



range of y 


99 -> 99 


gq -> gq + gq -> #g 


TOTAL 


all y 


5968.9 


1558.7 


7527.6 


\y\ < 0.5 


504.9 


129.3 


634.2 



Table 22: p-g (E T ) [mb GeV] 



range of 7/ 


99 -> 


gq ^ gq + gq ^ gq 


TOTAL 




2/ 


< 0.5 


48.17 


12.22 


60.39 



With the weaker gluon shadowing one finds 



Table 23: J 


dW/dy /or ^ = 5.5 ATeV 


range of y 


gg -> 


gq ^ gq + gq gq 


TOTAL 


all y 


24919.1 


3867.8 


28786.9 


\y\ <0.5 


2643.2 


438.4 


3081.6 



Table 24: gg (£ r ) [mb GeV] 



range of 7/ 


^ -» gg 


gq ^ gq + gq ^ gq 


TOTAL 




2/ 


< 0.5 


245.92 


40.95 


286.87 



A direct comparison between the results for shadowed and unshadowed 
parton distribution functions is shown in figure |TTJ and the pt distributions 
for LHC are shown in figure IT2. 
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Figure 11: Comparison of the rapidity distributions with strong (left figure) 
and weak gluon shadowing (right figure) to unshadowed distribution for LHC. 
The solid lines give the total contribution, the dotted ones depict the gg — > gg 
process and the dashed ones stand for the gq — > gq+gq — > gq processes. The 
thick lines again give the unshadowed results. 




23456789 10 123456789 

Pt Pt 

Figure 12: pt distributions for the two shadowing parametrizations at midra- 
pidity for LHC. 
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Therefore we find the following numbers at LHC: for unshadowed parton 
distributions one has at midrapidity 4845 gluons that carry a transverse 
energy of 16.4 TeV. The number density thus is n g = 363 fm 3 and the 
energy density is given by e s = 1229.7 GeV/fm 3 . The initial tempera- 
ture of an ideal gas derived with these numbers is Ti = 1169 MeV and 
T/^ =L0 = 1056.5 MeV for K=l. With the strong gluon shadowing we find 
634 gluons carrying a transverse energy of 1.93 TeV. We therefore have 
n g = 47.5 fm 3 and e s = 144.8 GeV/fm" 3 resulting in T ; = 684.9 MeV 
for K=1.5 and Tf= 10 = 618.9 MeV. With the shadowing version of |L5| 
we find 3082 gluons which carry a total transverse energy of 9.18 TeV, 
n g = 230.9 fm 3 , and e s = 678.6 GeV/fm 3 which results in a tempera- 
ture T; = 1011.08 MeV for K=1.5 and Tf= 10 = 913.62 MeV for K=l. 

5. Entropy production and ti multiplicities 
As is known, total entropy and entropy density, respectively, play a very im- 
portant role in the formation of a quark-gluon plasma. Total entropy reaches 
its final value when the system equilibrates and can, if assuming an adiabati- 
cal further evolution, be related to the effective number of degrees of freedom 



in the quark-gluon and in a pure pion plasma via |20], |2l 

s w (T c ) 



sw(T c ) 



0.7 ±0.2 (18) 



where s w and s qa are the entropy densities in the pion and quark-gluon 
plasma. The total entropy can then be related to the pion multiplicity as 

dy V d V J b =o r \ dy J b=0 

where c q9 = 4.02 for N f = 4 and c n w 3.6. 

A note on the separation between hard and soft processes is appropriate at 
this point. As emphasized above we introduced a cut-off at po = 2 GeV to 
ensure the applicability of perturbative QCD. Nevertheless there is always a 
soft component contributing to the production of transverse energy neglected 
in our studies so far. In it was shown that with p = 2 GeV at SPS the 
hard partons only carry about 4% of the total transverse energy Et- At RHIC 
energies they carry w 50% and for y/s = 2 TeV the hard partons already 
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carry « 80% of the total transverse energy. Since we here solely want to 
investigate the role of shadowing in hard reactions we will not calculate the 
pion multiplicity for RHIC where the soft contribution still is significant but 
restrict ourselves to the pion number at y ~ for LHC energies. 
If we employ the numbers for the entropy densities in the different plasmas 
and use our findings on the contributions of shadowing to the number of 
minijets we find that at y — one has 



, d y , 
, d y , 



6=0 



6=0 




3786, 

2413, (20) 
« 330, 



6=0 



when employing no shadowing, the '98 version of Eskola et al, and the strong 
gluon shadowing parametrization. 

6. Conclusions 

In this paper we investigated the influence of nuclear shadowing on ra- 
pidity spectra, transverse energy production and on macroscopic quanti- 
ties such as the initial temperature. We employed two different versions of 
parametrizations for the shadowing: one with a strong initial gluon shad- 
owing and a model independent one recently published by Eskola et al 



HJ. We found that the latter one gives an enhancement of minijet pro- 
duction at RHIC in contrast to the other case were a reduction to ~ 65% 
results. This difference directly manifests itself in the initial temperature 
Tj which happens to be smaller only for the strong gluon shadowing. At 
LHC the situation changes since there also the weakly shadowed gluons fi- 
nally result in lower spectra and Tj. Since the two shadowing parametriza- 
tions differ so drastically one finds a large difference in the results for the 
number of minijets at midrapidity: for the strong shadowing one has ~ 
630 gluons whereas for the weaker shadowing one finds ~ 3000 gluons. 
Since there are so few gluons for the strong gluon shadowing we find that 
the initial temperature at LHC is not dramatically higher than at RHIC! 
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